1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.00 1.01
## LDevsum 1.06 1.17
## dh0 1.72 2.67
## dh1 1.88 2.99
## dl0 4.07 8.16
## dl1 2.50 4.81
## dl2 1.60 2.48
## dl3 2.88 5.05
##
## Multivariate psrf
##
## 3.64
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1403.362 | 23642.734 |
| DIC3 | 1276.180 | 20946.750 |
| PWAIC | 103.750 | 1052.198 |
| WAIC | 1355.716 | 21358.294 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 100,000 iteration;
- 90,000 burn-in;
- 10 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.00 1.02
## LDevsum 1.00 1.01
## dh0 1.20 1.59
## dh1 1.19 1.58
## dl0 1.28 1.74
## dl1 2.04 3.67
## dl2 1.16 1.42
## dl3 2.43 5.08
##
## Multivariate psrf
##
## 2.1
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1448.01985 | 23921.990 |
| DIC3 | 1292.89309 | 21168.527 |
| PWAIC | 94.59621 | 1536.806 |
| WAIC | 1360.96116 | 21915.208 |